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Abstract 

Transition-metal-oxide based resistance random access memory is a promising candidate for 
next-generation universal non-volatile memories. Searching and designing appropriate new mate¬ 
rials used in the memories becomes an urgent task. Here, a new structure with the Ta02 formula 
was predicted using evolutionary algorithms in combination with first-principles calculations. This 
new structure having a triclinic symmetry (T'-Ta 02 ) is both energetically and dynamically more 
favorable than the commonly believed rutile structure (i?-Ta02). Our hybrid functional calcula¬ 
tions show that r-Ta02 is a semiconductor with a band gap of 1.0 eV, while i?-Ta02 is a metallic 
conductor. This large difference in electrical property makes Ta02 a potential candidate for resis¬ 
tance random access memory (RRAM). Furthermore, we have shown that T'-Ta02 is actually a 
Peierls distorted R-Ta02 phase and the transition between these two structures is via a directional 
displacement of Ta atoms. The energy barrier for the reversible phase transition from R-Ta02 
to r-Ta02 is 0.19 eV/atom and the other way around is 0.23 eV/atom, suggesting low power 
consumption for the resistance switch. The present findings provide a new mechanism for the 
resistance switch and will also stimulate experimental work to fabricate tantalum oxides based 
RRAM. 
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Faster and denser memories are ever-increasing demands and reqnirements of the ongoing 


advances in electronic devices and mnltimedia applications, etc . One of the most promis¬ 
ing candidates for next generation non-volatile memories is resistance random access mem¬ 
ories (RRAM), where data storage is achieved by the switch between a high resistance state 
(logic 0) and a low resistance state (logic 1) in certain materials by applying electric pnlses. 
The advantages of RRAM over the cnrrent commercial silicon-based Flash memory and other 
memory types are its simple strnctnre, good scalability, fast switching speed and good com¬ 
patibility with the widely nsed complementary metal-oxide-semicondnctor (CMOS) tech¬ 


nologies [2j . In a RRAM device, the choice of recording materials is important to achieve the 
above performances, and hence, searching and designing novel materials snitable for RRAM 
are received considerable focns. Nnmerons materials have been demonstrated to have the 
resistance switch behavior, among which the most promising type is transition-metal-oxides 
(TMO), snch as HfOa, j^, TiOj,^], TaOa, js-?] and so on. The typical RRAM device in¬ 
volving TMO has a sandwich-like Metal/TMO/Metal strnctnre, where the TMO layer is not 
restricted to one single type of oxide, for instance, the Pt/Si02/Ta205_a:/Ta02-x/Pt device 
has been recently developed by Samsnng Electronics j^. However, the microscopic mecha¬ 
nism of resistance switch in TMO-RRAM devices is very complex and not easy to nnravel. 


Previously, various mechanisms have been proposed [l|, J including filamentary conduction, 
space charge and traps, charge transfer, ionic conduction, etc. The filamentary conduction 
mechanism induced by the formation/migration of point defects has been reported in many 


TMO-RRAM 


a 


in| . For example, the diffusion of oxygen vacancies was found responsible 


for the resistance switch behavior in Ta205_a:/Ta02-x, where the resistance of the intrinsi¬ 
cally insulating Ta205 can be tuned by the generation of oxygen vacancies and clustering of 
Ta ions and while the metallic layer Ta02-a; mainly serves as the oxygen vacancy/ion reser¬ 
voir jol. Nevertheless, the complex microscopic mechanism of resistance switch in RRAM 
limits its practical applications. For large scale applications, obviously, a simple and con¬ 
trollable microscopic resistance switch in recording materials is vital. This is also the dream 
of materials scientists and a hot topic in RRAM research field. 


The above mentioned Ta02 was treated as a metallic material. However, in this work, 
we report a new stable triclinic structured Ta02 (r-Ta02) semiconductor with a band gap 
as large as 1.0 eV, in contrast to the metallic nature of the Ta02 with rutile structure that 
people normally believed. This new phase can be reversibly switched with its polymorph 
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conducting phase (i?-Ta02) via Peierls distortion where the phase transition energy barrier 
is 0.19 eV/atom from i?-Ta02 to r-Ta02 and is 0.23 eV/atom for the other way around. 
This huge resistance difference and easy phase transition between the two Ta02 polymorphs 
make Ta02 a promising candidate for RRAM. Furthermore, the present results provide a 
new resistance switching mechanism for RRAM, i.e., the metal-insulator phase transition 
via Peierls distortion, thus this work also expands the catalog of RRAM materials. The 
present hndings open a new way to design novel Ta02-based RRAM devices with lower 
power consumption, where the microscopic resistance switching is simple and controllable. 

Results 

Geometry and properties of the newly found T-Ta02. Since some substoichiomet- 
ric oxides exist or might be formed in the resistance-switch process of RRAM, the structure 
and property of any crystalline or even amorphous compounds formed by the transition 
metal and oxygen are vital for investigating the switch mechanism and designing novel 
RRAM devices. Here, using tantalum oxide as an example, we make a thorough exploration 


of the structures of oxic 


structure (R-Ta02) 11 


e Ta(IV)02. Ta02 is generally believed to have a rutile (tetragonal) 


13 l |. which is more stable than another symmetry I 4 i/a l 4 ]- The 


present evolutionary structure search and ab initio calculations (computational settings can 
be found in the Method section) revealed a new stable triclinic Ta02 phase (symmetry group 
F-i), for which the lattice parameters are: a = 90 . 0 °, /5 = 87 . 3 °, 7 = 95 . 7 °; a = 4.70 A, b 
= 5.20 A, c = 5.72 A. Details of the atomic structures can be found in the supplementary 
material. This newly found geometry contains 12 atoms in one unit cell, i.e., 4 Ta02 formula 
units, as shown in Fig. [H To make a clear comparison, a 1 x 1 x 2 supercell consisting of 12 - 
atoms for i?-Ta02 (symmetry group P 42 /mnm) was also shown. The lattice parameters for 
this i?-Ta02 supercell are: a = = 7 = 90 °; a = b = 4.99 A, c = 5.76 A. It is obvious that 

the most distinct difference between these two structures is the periodicity in the c direction, 
where for r-Ta02 it is almost twice that of i?-Ta02. While the deviation in lattice parame¬ 
ters along the a and b directions is within 6%. Further close analysis of the atomic positions 
in the two structures shown in Fig. [T] reveals that, besides the distortion of the octahedron 
formed by six O atoms, the labeled Ta atoms have the most signihcant displacements and 
these two inequivalent atoms displace almost half length of lattice parameter a. This large 
displacement results in a doubled periodicity along the c direction for r-Ta02. Therefore, 
r-Ta02 can be regarded as a distorted structure from the high-symmetric i?-Ta02. The 
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FIG. 1. (color online) Left: Unit cell of r-Ta02; Right: 1x1x2 supercell of R-Ta02. Blue and 
red balls represent Ta and O, respectively. Ta atoms showing the most significant displacements 
in the two structures are labeled as Tal and Ta2. 

phase transition from i?-Ta02 to r-Ta02 is induced by a typical Peierls distortion. 

The phase stability of the two Ta02 polymorphs are studied by the calculated cohesive 
energy and their phonon dispersions. Based on our PAW-GGA calculations, the cohesive en¬ 
ergy of r-Ta02 is 0.12 eV per formula unit lower than i?-Ta02. While our hybrid functional 
calculations showed that the cohesive energy of r-Ta02 is 0.20 eV per formula unit lower 
than i?-Ta02. This indicates that the exchange effect of electron interactions is quite differ¬ 
ent for these two phases with the same composition. Nevertheless, r-Ta02 is energetically 
more favorable than i?-Ta02. The results also clearly show that the structure distortions 
of i?-Ta02 lead to a more stable geometry. The phonon dispersions shown in Fig. [ 2 ] clearly 
unravel that r-Ta02 is also dynamically stable, while i?-Ta02 is not as shown by the imag¬ 
inary frequencies spreading over a large range of reciprocal space in Fig. [ 2 ](b). Therefore, 
the present found r-Ta02 phase is both energetically and dynamically more stable than 
i?-Ta02. 

To investigate the electronic property of the two Ta02 polymorphys, the partial density 
of states (PDOS) were calculated using hybrid functional, the results of which are shown 
in Fig. [31 It is clear that r-Ta02 is a semiconductor with a band gap as large as 1.0 
eV, while i?-Ta02 is metallic with no band gap. Furthermore, for both structures, the top 
of the valence-band mainly consists of the overlapping states of Ta-d and 0 -p, where the 
difference is whether these overlapped states are localized (for r-Ta02) or delocalized (for 
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FIG. 2. Phonon dispersion for (a) T-Ta02 and (b)ii-Ta02, the black dashed line indicating the 
position of zero frequency. 
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i?-Ta02). This significant difference in the electronic property of these two polymorphs 
makes Ta02 a potential candidate material for RRAM devices, where r-Ta02 is a high 
resistance state (HRS) and i?-Ta02 is a low resistance state (LRS). As mentioned above, 
T-Ta02 can be regarded as a distorted R-Ta02. This distortion leads to a more stable and 
donbled periodicity in the c direction and hence results in the metal-insulator transition 
(MIT), which we believe is a typical Peierls transition, similar to our previous work in GeTe 
system 

Phase transition between i2-Ta02 and T-Ta02 and the induced switching 
process. For an efficient RRAM device, as required by the read/write speed and power 
consumption, switch between high and low resistance states should be fast and easy. To 
explore the possibility of r-Ta02 to i?-Ta02 transition being the new candidate for a RRAM 
device, we have calculated the phase transition process, and the phase transition path and 
energy barrier are illustrated in Fig. HI Note that the hrst intermediate state next to R- 
Ta02 was found to be slightly more stable than i?-Ta02, showing distortion induced stability 
which is consistent with its dynamically unstable nature. Interestingly, as shown in Fig. HI 
the transition path between i?-Ta02 and r-Ta02 has two distinct periods, separated by a 
meta-stable state with its structure also shown in the hgure. Compared the metastable phase 
with the two end structures shown in Fig. [H especially the positions of the two numbered 
Ta atoms, we arrive at the conclusion that these two Ta atoms play the pivotal role in the 
MIT where they move independently but directionally with the largest migration distance 
to accomplish the phase transition between i?-Ta02 and r-Ta02. During the transition 
from R-Ta02 (LRS) to r-Ta02 (HRS), atom Ta 2 hrstly migrates a large distance to the 
similar symmetry position as in r-Ta02 leading to a meta-stable state (the inset of Fig. H]), 
where the energy barrier to reach this state is 0.14 eV/atom. Then it is followed by a large 
migration of atom Tal leading to the structure of r-Ta02, where the energy barrier for this 
meta-stable phase transforming to r-Ta02 is 0.11 eV/atom. The overall barrier to overcome 
for the MIT transition from R-Ta02 to r-Ta02 is around 0.19 eV/atom. The other way 
around is for the phase transition from r-Ta02 (HRS) to i?-Ta02 (LRS) as applying an 
opposite electric held, where the overall energy barrier is 0.23 eV/atom. At this stage, it is 
obvious that Ta02 is a suitable candidate material for RRAM devices, where the reversible 
switch between HRS and LRS is a Peierls distortion type which can be achieved by applying 
different electric helds of opposite directions and various energies. We believe this new Ta02- 
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FIG. 3. Partial density of states for (a) r-Ta02'^nd (b)i?-Ta02. The vertical dashed lines indicate 
the position of Fermi level. 










































RRAM should be low power consumption. Even though there is no direct relation between 
the energy barrier of this metal-semiconductor transition with the power consumption of real 
devices, it is reasonable that the larger the energy barrier, the higher the power consumption. 
The conclusion is obvious as compared Ta02 with other materials also used in random access 
memory. For example, in the well-known Ge-Sb-Te (GST) system for phase-change random 
access memory (PGRAM) where the reversible phase transition between amorphous (HRS) 
and cr ysta lline (LRS) under various electric pulse is used for data storage 


16 


et ah 


191 , Kolobov 


20 | recently reported that the energy barrier between the ordered-crystalline and 


the distorted structure leading to the amorphous state is as large as 0.45 eV/atom, almost 
twice that of the barrier in our study. Therefore, we believe that the power consumption 
when using the metal-semiconductor transition between the two Ta02 polymorphs in RRAM 
could be lower than that of GST in PGRAM. Moreover, like the hlamentary mechanism that 
commonly exists in the RRAM based on TMO, the resistance switch between i?-Ta02 and 
r-Ta02 normally occurs in a tiny area of the material in practical RRAM devices, not 
necessarily inside the whole bulk. This will result in an even lower energy consumption and 
faster data-storage speed, i.e., make the device ’greener’. Finally, it is worth to mention that 
the above calculated negative phonon frequencies in i?-Ta02 will not affect the application 
of this material in RRAM. A similar case is in Ge2Sb2Te5, a prototype material for PGRAM, 
where the resistance switch is accomplished by the reversible transition between a metastable 
rocksalt structure and an amorphous state. For metastable rocksalt structured Ge2Sb2Te5 
negative phonon frequencies have also been reported by theoretical calculations 21|. 

Conclusion 

In summary, by using evolutionary algorithms, we found a new tantalum dioxide {T- 
Ta02) structure, which is more stable than the commonly believed rutile symmetry {R- 
Ta02). The newly found r-Ta02 which can be regarded as distorted R-Ta02 is a semicon¬ 
ductor, in contrast to the metallic nature of i?-Ta02. This Peierls distortion induced re¬ 
versible metal-semiconductor transition with low energy barriers of 0.19 and 0.23 eV/atom, 
respectively, provides a new mechanism for RRAM devices and expands the catalog of 
RRAM materials. Finally, our work will stimulate the re-building of transition-metal-oxide 
phase diagram aiming at Ending new structures with various electronic properties, and also 
promote the experimental work to fabricate RRAM devices based on these transition metal 
oxides. 
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FIG. 4. Transition path from i?-Ta02 to r-Ta02, the energy difference between various states 
is labeled in the unit of eV/atom. Structure of the meta-stable intermediate state as denoted by 
the solid arrow is also shown, where the Ta atoms corresponding to numbered ones in Fig. [1] are 
labeled as Tal and Ta2. 


Methods 

Structure prediction. Evolutionary algorithm as implemented in the code USPEX 
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231 ] was used to predict the structures with the formula Ta02. During the calculations, the 


already known rutile Ta02 structure was introduced as the seed. 

Ab initio calculations. All the structures screened by USPEX were fully relaxed using 


VASP code 


24 


26 | . in conjunction with projector augmented wave potentials within the 


generalized-gradient approximation (PAW-GGA) ( 27 |, while for the calculation of the band 
gap of structures, hybrid functional HSE 06 was employed to deal with the localized d 
electrons of Ta. Pseudopotentials with electronic conhgurations of and 2 s^ 2 p‘^ were 
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used for Ta and O, respectively. The structures mentioned in this paper were fully relaxed 
until the forces on them were less than 10 meV/A with a cutoff energy 450 eV. i^-point 
meshes for structures with different symmetry and lattice size were generated with the same 
resolution. Our calculations also showed that the spin-polarization effects in the structures 
involved in the present study were negligible. Phonon dispersions were calculated by Density 
Functional Perturbation Theory (DFPT) as implemented in the code Phonopy 29 |]. The 
phase transition process was investigated using the generalized solid-state Nudged Elastic 
Band (G-SSNEB) method 30 |. First, we assumed a linear transition path by sampling 8 
intermediate states between the two end-structures (i?-Ta02 and r-Ta02), and during the 
Minimum Energy Path (MEP) optimization, besides the atomic positions, the shape of the 
supercells is also allowed to relax. 
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